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Abstract 

The propagation of the isoscalar and isovector sound modes in a hot nuclear 
matter is considered. The approach is based on the colhsional kinetic theory 
and takes into account the temperature and memory effects. It is shown that 
the sound velocity and the attenuation coefficient are significantly influenced 
by the Fermi surface distortion (FSD). The corresponding influence is much 
stronger for the isoscalar mode than for the isovector one. The memory 
effects cause a non-monotonous behavior of the attenuation coefficient as a 
function of the relaxation time leading to a zero-to-first sound transition with 
increasing temperature. The mixing of both the isoscalar and the isovector 
sound modes in an asymmetric nuclear matter is evaluated. The condition 
for the bulk instability and the instability growth rate in the presence of the 
memory effects is studied. It is shown that both the FSD and the relaxation 
processes lead to a shift of the maximum of the instability growth rate to the 
longer wave length region. 

PACS number: 21.60.Ev, 21.65. -hf, 24.10.Nz 
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I. INTRODUCTION 



In the vicinity of the equihbrium state nuclear matter is stable with respect to parti- 
cle density and Fermi surface distortions and the excitation of both the isoscalar and the 
isovector sound modes is possible. Propagation of an isoscalar sound wave in nuclear matter 
depends crucially on the Landau parameter Fq in the quasiparticle interaction amplitude 
and on the relaxation processes. For zero temperature, there is an underdamped zero-sound 
mode at Fq > because the phase velocity of the sound wave exceeds the speed of particles 
inside the Fermi sphere. A strong Landau damping appears at — 1 < Fq < where a nonzero 
transfer of energy from the wave to particles is possible |^ and the wave transforms to an 
overdamped mode. The picture of propagation of the zero-sound wave is essentially more 
complicated in the case of a hot nuclear matter. Due to the existence of the temperature 
tail of the equilibrium distribution function, the phase conditions for the Landau damping 
are fulfilled here for positive values of the Fq and a possibility for a propagation of the sound 
wave appears in the region — 1 < Fq < 0. 

The zero sound is transformed to the first-sound mode in the limit of strong interaction, 
|Fo| S> 1, [0] or in the frequent collision regime at high temperatures |Q. It is necessary to 
stress that the sound velocity c is directly related to the nuclear matter incompressibility 
coefficient K for the first-sound mode only. In this case, one has c ~ Ci = ^K/9m, where ci 
is the velocity of the first sound. In general, the sound velocity c is a complicated function 
of both K and the dimensionless coUisional parameter ut, where r is the relaxation time 
and u is the eigenfrequency of the sound mode. In the present work we obtain a simple 
analytical expression for the sound velocity c which provides a description for both the 
frequent- and rare-collision limit as well as for the intermediate cases. Special attention is 
paid to the propagation of the isovector sound and for the zero- to first- sound transition 
for this mode. It is well known that hydrodynamic approaches like the Goldhaber- Teller 
or Steinwedel- Jensen models give a reasonable description of the nuclear isovector giant 
resonances (IVGR). However it is not the case for the isoscalar giant resonances (ISGR), 
where the Fermi surface distortion effects plays an important role 1^-0 and the traditional 
hydrodynamic model can not be applied. In our approach, this situation is directly related 
to the peculiarities of the propagation of both isovector and isoscalar sounds. 

The interparticle collisions on the distorted Fermi surface lead to the coUisional damping 
of the sound wave. Two limiting regimes ut and ut oo provide the existence of 
the non-damped first and zero sounds, respectively, lU. In what follows we will use the 
coUisional kinetic theory, taking into account the memory effects on the collision integral 
p|JT0HT2| . We will discuss a special feature of the temperature dependence of the attenuation 



of the sound mode in hot nuclear matter. In particular, we will show that memory effects 
lead to a bell-shaped form of the attenuation coefficient k as a function of the temperature, 
providing a correct behavior of k in both limiting regimes ujt —>■ and ut — *• oo. 

With decreasing bulk density or increasing temperature the nuclear matter reaches the 
regions of mechanical or thermodynamical instabilities with respect to small particle density 
distortions P,p!3HT7| and to separation into liquid and gas phases [ p!8| , p!9[ . The general 



instability condition of the Fermi liquid reads 1 + Fk/{2k -|- 1) < [|20|, where Fk is the 



Landau's parameter in the expansion of the quasiparticle interaction amplitude in Legendre 
polynomial pi] . However the development of instability depends not only on the equation of 
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state, but also on the dynamical effects such as the dynamical Fermi-surface distortion (FSD) 
effect [^. The FSD effects strongly reduce, by the factor ~ (|Fo| — 1)^''^, the instability 



growth rate F in the unstable region < — 1 — Fq -C 1 0. In the present work, we will 
consider the influence of the Fermi-surface distortion, relaxation and memory effects on the 
instability growth rate F. 

The plan of the paper is as follows. In Sec. II we derive the general equation of motion 
for the particle density vibrations in the presence of the memory effects. We introduce the 
renormalized incompressibility coefficient and the viscosity coefficient which are frequency 
dependent due to the memory effect. In Sec. Ill we obtain a simple expression for the 
refraction and attenuation coefficients for both the isoscalar and the isovector sound waves 
which propagate in a hot nuclear matter. We consider also the mixing of the isoscalar and 
isovector modes in an asymmetric nuclear matter taking into account the relaxation and the 
memory effects. The development of the bulk instability in a low density nuclear matter 
in the presence of the relaxation and the memory effects is considered in Sec. IV. The 
conclusion is given in Sec. V. 



II. SOUND PROPAGATION IN HOT NUCLEAR MATTER 

We start from the collisional kinetic equation 

df pdf dVdf 

^ = (1) 

ot m or or Op 

Here (5St(t) is the collisional integral, V = V{r,t) is the self-consistent mean field and / 
= /(r, p, t) is the Wigner distribution function in which we will take into account only the 
distortion of the Fermi sphere with multipolarities i < 2 

f = fs + Sf, 6f = j26U (2) 

i=i 

The distribution fs = /s(r, p, t) corresponds to the spherical Fermi surface and 6f represents 
both the quadrupole deformation and the displacement of the Fermi surface. For small 
deviations from a Fermi sphere the right-hand side (RHS) 5St(i(:) of ([I|) is a collision integral 
linearized in Sf and it can be represented in the form 

(5St(t) = r dt'A{t-t')6f{t'), (3) 

J —oo 

which takes into account the memory effects due to the memory kernel A(t — t'). In this 
paper we will not use an explicit form of A. Below we will need (5St(t) only for periodic 
oscillation of 6f with the eigenfrequency u. Assuming the restriction i < 2 for the Fermi 
surface distortion, the collision integral (|^) can be written in the form of the extended 
r-approximation, see Refs. P^P3| , 

SSt{t) = (4) 
where the relaxation time t^^ is cu-dependent due to the memory effects. 
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The assumption in Eq. allows us to reduce the collisional kinetic equation (0) to 
the local equations of motion for the particle density p = p{r,t), the displacement field 
Xa = Xali")^) and the pressure tensor Pa/3 = Pa/3{r,t) derived as the lowest p-moments of 
the distribution function /(r, p, t): 

9 dp dxa I f g dp Pa 1 r g dp 



{27:%^ dt pJ '"^ m J {2T:hf 

(5) 

Here Ua = dxajdt is the velocity field and g is the spin-isospin degeneracy factor. Taking 
the first two p-moments of Eq. (|l]) one obtains, see Refs. p^-p4|, 



where P is the pressure due to motion of nucleons without distortion of the Fermi sphere and 
Pl^ is associated with quadrupole distortion of the Fermi surface: Pq/? = P(5q/3 + P^^. The 
pressure tensor Pl^ is responsible for the dissipative processes. Taking the second p-moment 
of Eq. (|l|) one obtains the following equation for the pressure tensor 

Q Q d d 

The equations of motion (H) and (|^ are closed. They can be applied to both the isoscalar 
and isovector sound excitations. Let us consider the isoscalar compression mode. To simplify 
the problem, we can rewrite the expression in the parenthesis in Eq. (P) near the equilibrium 
value of the density pcq as 



9 p d ^ d 6e d 



V 



bp 



(8) 



bp' 

where index "eq" refers the equilibrium state, e is the energy density of particles 

e = Ckin + Cpot, (9) 

ekin is the kinetic energy density 

''-=2'^ = Jj2^^2^^^ = -0m^T^ P 

and epot is the potential energy density which is related to the mean field V by 

V = be^ot/bp. (11) 

Note that {be/bp)eq is the chemical potential, which does not depend on the space coordinate 
r, for the equilibrium state of the nucleus. We have used this fact when deducing Eq. 
(I). Solving Eq. (0) with respect to P^^, using Eq. (||) and the continuity equation 
bp = p — Pcq = — div(peqX)) 'WG obtaiu an equation for the density vibration in the form 
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u'^dp + {KjQm)V'^5p = icj(4r/^/3mpeq)V^(5p. (12) 

Here 

K = K + 8(e,,„/p),Jm f T^f^) , (13) 
where K = 9(5^e/5p^) peq is the static incompressibihty and r]^ is the viscosity coefficient 



^^ = Re( ^ IPeq. (14) 

We point out that there is a significant difference between the static nuclear incompressibihty 
coefficient, K, i.e., derived as a stiffness coefficient with respect to a change in the bulk 
density, and the dynamic one, K'^ of Eq. (0), associated with the sound propagation. This 
difference is due to the second term on the RHS of Eq. ([T3|) caused by the Fermi-surface 
distortion effects. The quantity rj^^ in Eq. (|T^ determines the time irreversible contribution 
to the pressure tensor P^^ and can be considered as the viscosity coefficient due to the 
relaxation occurring on the distorted Fermi surface. Expression (113) is valid independently 



of the nucleon's collision rate. The viscosity goes to zero in both the rare, Tr^^ — > oo, and 
frequent, r^^^ 0, collision limits. 



III. DISPERSION RELATION, MEMORY EFFECTS AND DAMPING 



Assuming a plane wave solution Sp ~ exp(iq- r—iut) one obtains from Eq. (12) the 
following dispersion relation 

^2 = «/9m)g^ - z^(4r/^/3mpeq)gl (15) 

The solution of this equation defines the complex wave number q [u is real). A simple 
solution to Eq. (|l^) can be obtained in two limiting cases of the frequent collision (first 
sound) regime, uTr^^ — > 0, and the rare collision (zero sound) regime, uTr^uj — ^ oo. The sound 
velocity c = u/q is given by 



Cl 



K/9m if LUTr^uj and c = Cq = \J {K + Afr)/9m if uJTr^i^ oo, (16) 



where AK fa 8(ekin/p)eq ~ (24/5) e^? ~ 200 MeV (we adopted the kinetic Fermi energy ep ^ 
40 MeV). We point out that the value of AK is comparable with the static incompressibihty 
K pa 220 MeV and we have cq ~ V^Ci. The factor a/2 in this relation is due to the restriction 
£ < 2 for the multipolarity £ of the Fermi surface distortion. In a general case of arbitrary 
i this factor is increased to -\/3 ||l]]. The result ([TBD means that in contrast to the first 
sound (frequent collision) regime, the sound velocity of the compression mode can not, in 
general, be used directly to extract the static incompressibihty of K because of the additional 
contribution from the Fermi surface distortion effects which result in the renormalization of 
the incompressibihty K into K'^. 

Using both asymptotic sound velocity ci and cq, the solution to the dispersion relation 
(|15|) can be written as 
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q = —{n + ik), (17) 
Co 

where the refraction coefficient n and the attenuation coefficient k (both real) are obtained 
from the following equation 



n + iK=. : . 18) 



In the frequent collision (first sound) regime we obtain from Eq. (|T8| 



n = — , K = UJTr^^{Co/2Ci)[{Co/Cif -1] if (19) 

Cl 

In the opposite case of the rare collision (zero sound) regime we obtain 

n = l, K=[l- {ci/Cof]/{2uJTr,^) if UJTr,^ > 1. (20) 

The attenuation coefficient k in both limiting regimes is a complicated function of the 
frequency u because of the memory effect in the relaxation time Tr^cu- In the case of sound 
propagation in hot nuclear matter, the competition between the temperature smoothing 
effects in the equilibrium distribution function and dynamical distortions of the particle 
momentum distribution leads to the following expression for the relaxation time [^|11],|12| 

where T is the temperature of nuclear matter and the ci;-dependence of r^^; is due to the 
memory effects in the collision integral. Below we will use ^ = I/Att"^ [|^ andro = ah, a = 
9.2 MeV (for the isoscalar mode) |26[] . 



Equations (^) and (|T8[) are valid for arbitrary collision times Tr^uj and thus describe both 
the zero and the first sound limit as well as the intermediate cases. From it one can obtain 
the leading order terms in the different limits mentioned. In Fig. 1 we have plotted both 
coefficients n and k as obtained from Eq. ([T8|). In the high temperature limit, the system 
goes to the frequent collision (first sound) regime with the saturated refraction coefficient n 
~ Co/ci ~ -\/3 (we use the factor y/S instead of \/2 assuming the contribution of the higher 
multipolarities i > 2 in the Fermi surface distortion as was mentioned above) and the 
attenuation coefficient k ~ Tr^uj ~ l/T^. In the opposite low temperature limit, the system 
is close to the zero sound regime with n ~ 1. We point out a shift of both n and n by 
nonzero values at T — 0. This is due to the memory effect in the relaxation time r^^i^ of Eq. 



(pip: in the very high frequency limit, the system can exist close to the first sound regime 
at n y/S even at zero temperature. The position of the maximum of k{T) in Fig. 1 can 
be interpreted as the transition temperature T^r of zero- to first- sound regimes in a hot 
Fermi system. The value of Ttr depends slightly on the sound frequency uj and it is shifted 
to smaller values with the increase of u. 

Let us consider now the isovector sound mode in a symmetric nuclear matter with Pn.eq = 
Pp^eq, where pn,eq and pp^eq are the equilibrium neutron and proton density respectively. The 
general equations of motion (|^) and (^ are still correct. However the energy density e in Eq. 
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is related now to the symmetry energy -E'symm- The corresponding first sound velocity ci 



for the isovector mode is given by, Ch. 6 



Ci = sj2Esjmm/m, (22) 

where i?symm = (1/3)6^7(1 + Fq) ^ 30 MeV, and Fq is the isovector Landau parameter in the 
quasiparticle interaction amplitude. The zero sound velocity cq for the isovector mode can 
be found from Eqs. (|]) and (|^) in the rare collision limit Tr^^j — ^ 00. Taking into account Eq. 
(p!T|), one obtains, see also Eq. ([T6|), 



Co = Y 2 (-E'symm + A-^symm) /"^, (23) 

where AEgymm ~ (4/9)(efci„/p)gq ~ (4/15) ep ^ 10 MeV. We point out that, in contrast to 
the isoscalar mode, the Fermi surface distortion effect leads to a relatively small increase of 
the isovector zero sound velocity cq, Eq. (|23|) , with respect to the first sound one ci, Eq. 
(p^). The dispersion relation (|I5|) takes the form 



= (2^symm,c^/"^)?^ " ^^(4r/^/3mpeq)g^ (24) 



where 



E'syram,. = ^symm + (4/9) (e,,„/p),qlm ^^.'''^ . (25) 

All relations (|l^)-(^) are still correct in the isovector case if both asymptotic velocity Cq 
and Ci are taken from Eqs. (p2|) and (^). In Fig. 2 we have plotted the coefficients n and k 
as obtained from Eq. ([T8|) for the isovector mode with the collision parameter a = 4.6 MeV 
26| . We point out that the transition temperature Ttr of zero- to first- sound regimes for 



the isovector mode is significantly smaller than Ttr for the isoscalar one. 

In an asymmetric nuclear matter, both the isovector and the isoscalar modes are depen- 
dent on each other. The particle density fiuctuation bp takes a bispinor form bp = {Sp+, 
Sp-), where Sp+ and 6p- are the isoscalar and isovector components respectively. A solu- 
tion of the corresponding equations of motion (|^) and (|^) leads to the following dispersion 
relation, see also Eqs. ([T5|) and (p4D, 



, - {K'^/9m)q^ + iiu{Ari^/3mp^^)q^ {Ix/m)q^ + 0{P , , _ 

' {Iy/m)q' + 0{P) - {Ei^^^j9m)q' + iu:{A^^/3mp,^)q' ) 

(26) 



Here, / = {pn — Pp)eq/ {pn + Pp)eci "C 1 is the asymmetry parameter, p„ and Pp are the neutron 
and proton densities respectively, peq = {pn + Pp)eq and the coupling constants x and y are 
given by 

x = -(2/9)(4e^-i^/3), = -(2/9)(4e^ + ir/6 - 9Esymm). (27) 

As is seen from Eq. (P?]), the eigenfrequency uj and the corresponding sound velocity for 
both the isoscalar and the isovector modes are independent of each other in the linear order 
of the asymmetry parameter /. The structure of bispinor 6p = {6p^, 5pJ) is different for 
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both the isoscalar-hke and the isovector-hke modes. For the isoscalar-hke mode with the 
eigenfrequency given by a solution to Eq. ([TsD , the main contribution to the bispinor 6p 
is due to the isoscalar component 5p+ ~ 1 while the isovector component (5p_ is proportional 
to the asymmetry parameter /. Namely, 



The opposite situation takes place for the isovector-like mode with the eigenfrequency tUp 
given by a solution to Eq. (123). In this case, one has 6p^ ~ 1 and 5p+ ~ /. 



IV. BULK INSTABILITY 

Let us consider now the bulk instability regime K < and introduce an instability 
growth rate F = — icj (F is real, F > 0), see Refs. [0,^. The amplitude of the density 
fluctuations, 6p ~ exp(iq ■ r—iujt) ~ exp(Ft), grows exponentially if F > 0. To prevent an 
unphysical infinite growth of the short wave length fluctuations (see dotted line in Fig. 3), 
we will take into account the velocity dependent contribution to the effective interparticle 
interaction. Due to the corresponding change in the selfconsistent mean field V in Eq. (^, 
an additional anomalous term ~ appears in the dispersion relation ([T5|) and the equation 



for the instability growth rate F takes the following form |2T[ 

T' = \c,\'q^-aT)q'-Ksq\ (29) 



Here ci = iJ\K\/Qm and 



C(r) = TT^T^e^, (30) 
3 m 1 + Fr.. 



where is the relaxation time = ah/T"^ [Q. The constant Kg in the anomalous dispersion 
term in Eq. ( PPD depends on the model. In the case of the effective Skyrme forces, one has 
pT| Kg = fi^/dm^ + {9ti — 5)f:2)po/32m, where ti and t2 are the parameters of the velocity 



dependent part of the Skyrme forces We point out that the instability regime with 



< can be reached at a low bulk density po- The incompressibility K is given by 

K = 6eF{l + Fo) {l + F,/3)-\ (31) 

Here, the Landau parameters are related to the parameters t„ of the effective Skyrme 
forces. Namely, 



9po 
Sep 





3 1 


m* 




fm* \ 




^ 2 ^3 Po 


— + 3 




, Fi = 3 1 
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V m J 


V m / 



(32) 



where m/m* = 1 + mpo (3ti + 5t2)/8h'^. For the commonly used set of parameters tn, the 
instability regime with Fq < —1 is reached at po ^ O.Speq. 

In Fig. 3 we have plotted the instability growth rate F as obtained from Eq. (pQ]). The 
calculation was performed for the Skyrme force SIII. For the relaxation time we used 
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a = 9.2 MeV and the bulk density po was po = 0.3 peq, where peq is the saturated density 
Peq = 0.1453 fm~^. We also show in Fig. 3 the result for the nonviscous nuclear matter 
neglecting the anomalous dispersion term ~ and the Fermi surface distortion effect (dotted 
line). The non-monotonous behavior of the instability growth rate as a function of the wave 
number q is due to the anomalous dispersion term in Eq. (|29|), induced by the velocity 
dependent terms in the interparticle interaction. The instability growth rate F reaches a 
maximum, Fmax, at a certain q = gmax and F goes to zero at g = gcrit- The existence of 
the critical wave number gcrit = Icip/^^ for an unstable mode is a feature of the system 
with the anomalous dispersion term |T3[. The distortion of the Fermi surface leads to a 



decrease of the critical value gcrit, i-e., the nuclear matter becomes more stable due to the 
FSD effect. We can also see that the presence of viscosity and the FSD effect lead to a 
shift of the position gmax of the maximum of F(g) to the left. Thus, the instability of the 
nuclear matter with respect to short-wave-length density fluctuations decreases due to the 
viscosity and the FSD effect and the most unstable mode is shifted to the region of the 
creation of larger clusters in the disintegration of nuclear matter. For a saturated nuclear 
liquid one has for the force parameters to < 0, ts > and ts > 0. However both values gmax 
and gcrit have a non-monotonous behavior as a function of the bulk density po because of 
the additional po— dependence of the Fermi energy ep in Eqs. ( PID and (^). The particle 
density dependence of the values of gmax and gcrit is shown in Fig. 4. The instability growth 
rate F(g) as well as the values of gmax and gcrit are only slightly sensitive to the change of 
temperature at T ^ 10 MeV, where the temperature dependence of the bulk density po can 
be neglected. A more sophisticated consideration is necessary near the critical temperature 
Tcrit ~ 17 MeV where the nuclear matter is unstable with respect to the liquid-gas phase 
transition. 



V. CONCLUSION 

Starting from the collisional kinetic equation (|l]), we have derived the dispersion relations 
(^) and (|29|) for both the stable and the unstable regime of the density fluctuations in a 
heated nuclear matter. The dispersion relations are influenced strongly by the FSD effect 
and the anomalous dispersion term. The presence of the Fermi surface distortion enhances 
the stiffness coefficient for a stable mode and reduces the instability growth rate for an 
unstable one. There is a significant difference between the static nuclear incompressibility 
coefficient, K, i.e., derived as a stiffness coefficient with respect to a change in the bulk 
density, and the dynamic one, K'^ associated with the zero sound velocity, see Eq. (^3[). 
The FSD effect is responsible for the collisional relaxation of the collective modes in the 
Fermi liquid and for the non-Markovian character of the nuclear matter viscosity (memory 
effect in the viscosity rj^, Eq. (|T4|)). The memory effects in the viscosity play an essential 
role in the description of the temperature dependence of the refraction coefficient n, see Eqs. 
(p!8|) and (|2l|) . We have noted also the bell-shaped form of the attenuation coefficient k as a 
function of the temperature T, see Figs. 1 and 2 . This pecuharity of k(T) provides a new 
criterion for the determination of the transition temperature Ttr between the zero-sound and 
first-sound regimes in hot nuclear matter. 

Our consideration provides a good basis for understanding the difference in the devel- 
opment of the spinodal instability in nuclear matter taking into account both the FSD and 
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the viscosity effect. We liave sliown tliat tfie instability growtli rate T{q) in an unstable nu- 
clear matter with velocity dependent effective interparticle interaction is a non-monotonous 
function of the wave number q because of the anomalous dispersion term. The anomalous 
dispersion term removes an unphysical infinite growth of the short wave length fluctuations. 
The non-monotonous behavior of the instability growth rate T{q) is characterized by two 
wave numbers gmax and gcrit- We point out that both the FSD effect and the relaxation 
processes lead to a shift of gmax to the longer wave length region, providing an increase of 
the relative yield of heavier clusters. 

The main results were obtained assuming a quadrupole distortion of the Fermi surface. 
We point out, however, that the key expression (0) for the viscosity coefficient can be 
rewritten identically in the following form 



where ci and cq are the first and zero sound velocity respectively. This expression can also 
be established from a general consideration, see Ref. [|29| Ch.8, and can be used for arbitrary 
multipolarities of the Fermi surface distortion. 



This work was supported in part by the US Department of Energy under grant ^ FG03- 
93ER40773. We are grateful for this financial support. One of us (V.M.K.) thank the 
Cyclotron Institute at Texas A&M University for the kind hospitality. 




(33) 



VI. ACKNOWLEDGEMENTS 



10 



REFERENCES 



[1] E. M. Lifshits and L. P. Pitaevsky, Physical kinetics, (Pergamon Press, Oxford - New 
York - Seoul - Tokyo, 1993). 

[2] V. M. Kolomietz, A. B. Larionov and M. Di Toro, Nucl. Phys. A613, 1 (1997). 

[3] C. J. Pethick and D. G. Ravenhall, Ann. of Phys. 183, 131 (1988). 

[4] A. A. Abrikosov and I. M. Khalatnikov, Rep. Prog. Phys. 22, 329 (1959). 

[5] J. M. Eisenberg and W. Grainer. Nuclear Models, Collective and Single-Particle Phe- 
nomena. North- Holland, Amsterdam, 1970, Ch.lO. 

[6] G. F. Bertsch, Nucl. Phys. A249, 253 (1975). 

[7] G. Holzwarth and G. Eckart, Nucl. Phys. A325, 1 (1979). 

[8] J. R. Nix and A. J. Sierk, Phys. Rev. C 21, 396 (1980). 

[9] V. M. Kolomietz, Sov. J. Nucl. Phys. 37, 325 (1983). 
[10] G. F. Bertsch, P. F. Bortignon and R. A. Broglia, Rev. Mod. Phys. 55, 287 (1983). 
[11] S. Ayik and D. Boilley, Phys. Lett. B276, 263 (1992); B284, 482(E) (1992). 
[12] V. M. Kolomietz, A. G. Magner and V. A. Plujko, Nucl. Phys. A545, 99c (1992). 
[13] C. J. Pethick and D. G. RavenhaU, Nucl. Phys. A471, 19c (1987). 
[14] M. Colonna, Ph. Chomaz and J. Randrup, Nucl. Phys. A567, 637 (1994). 
[15] M. Colonna and Ph. Chomaz, Phys. Rev. C 49, 1908 (1994). 
[16] S. Ayik, M. Colonna and Ph. Chomaz, Phys. Lett. B353, 417 (1995). 
[17] V. Baran, M. Colonna, M. Di Toro and A. B. Larionov, Nucl. Phys. A632, 287 (1998). 
[18] H. Jaqamann, A. Z. Mckjian and L. Zamick, Phys. Rev. C 27, 2782 (1983). 
[19] H. Miiller and B. D. Serot, Phys. Rev. C 52, 2072 (1995). 
[20] I. Ya. Pomeranchuk, Sov. Phys. JETP 8, 361 (1959). 
[21] V. M. Kolomietz and S. Shlomo, Phys. Rev. C 60, 044612 (1999). 
[22] V. M. Kolomietz, V. A. Plujko and S. Shlomo, Phys. Rev. C 52, 2480 (1995). 
[23] D. Kiderlen, V. M. Kolomietz and S. Shlomo, Nucl. Phys. Nucl. Phys. A608, 32 (1996). 
[24] V. M. Kolomietz and H. H. K. Tang, Phys. Scripta. 24, 321 (1981). 
[25] A. Kolomiets, V. M. Kolomietz and S. Shlomo, Phys. Rev. C 59, 3139 (1999). 
[26] V. M. Kolomietz, V. A. Plujko and S. Shlomo, Phys. Rev. C 54, 3014 (1996). 
[27] A. Bohr and B. R. Mottelson, Nuclear Structure, Vol. 2 (Benjamin, New York, 1975). 
[28] D. Vautherin and D. M. Brink, Phys. Rev. C 5, 626 (1973). 

[29] L. D. Landau and E. M. Lifshits, Fluid Mechanics, (Pergamon Press, Oxford, 1963). 



11 



FIGURES 




FIG. 1. Refraction, n, and attenuation, k, coefHcients of the isoscalar sound wave as functions 
of temperature. The calculation was performed for two eigenenergies Hu = IMeV (solid line) and 
hu = leV (dashed line). 
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FIG. 2. Same as Fig. 1 for isovector mode. 
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FIG. 3. Dependence of the instability growth rate T on the wave number q. The calculations 
were performed for Skyrme force SIII, temperature T = 6 McV and density po = x peq with x = 0.3 
and Peq = 0.1453 fm~^. The sohd curve is for the viscous nuclear matter with a = 9.2 MeV 
including both the memory and the Fermi-surface distortion effects. The dashed and dotted lines 
are the results for the nonviscous liquid: curve (1) is for a nuclear matter neglecting the FSD effect; 
curve (2) is the result in the presence of the FSD effect and dotted line is the same as curve (1) 
neglecting the anomalous dispersion term. 



14 



2 



(J 
D" 



X 
03 



15 - 



1 - 



05 



I 
I 
I 
I 

I / 



^max 







01 02 03 04 05 06 07 



FIG. 4. Dependence of the characteristic wave numbers of the instabihty growth rate r(g) on 
the dimensionless density parameter x = po/Peq- The calculation was performed for Skyrme force 
SIII and temperature T = 6 MeV. 
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